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1.  Introduction 

In  the  Arctic  Ocean  and  Its  surrounding  seas,  an  enormous  number  of  tee  floes  of  various 
sizes  and  shapes  are  floating  on  the  surface.  Most  of  these  ioe  floes,  under  the  Influence  of 
winds,  ocean  currents  and  tidal  waves,  are  In  constant  but  chaotic  motion  and  are  continuously 
being  deformed.  Any  such  floe,  individually  considered  over  a  discrete  period  of  time,  repre¬ 
sents  a  microscopic  description  of  ice  movement  in  the  Polar  Oceans.  On  the  other  hand,  for 
many  practical  reasons,  we  are  interestedtln  the  average  motion  of  ioe  over  a  large  area  such 
as  over  1000  kilometer  radius  or  larger,  and  for  an  extended  period  of  time,  say  several  days, 
months  or  years.  The  description  of  these  large  scale  movements  of  pack  ice  should  be  given 
by  the  average  of  local  ice  movement  over  a  large  area  and  an  extended  period  of  time.  This  is 
a  macroscopic  description  of  pack  ice  movement.  In  order  to  correlate  the  microscopic  and 
macroscopic  description  of  pack  ice  so  that  we  may  predict  the  pack  ice  deformation  and  move¬ 
ment,  the  development  of  a  mathematical  model  is  necessary.  There  are  several  mathematical 
models  in  literature  (Fel'zenbaum  1958,  Campbell,  1965,  Doronln  1970,  Coon  et  al,  1974  and 
Pai-Li  1977).  In  most  of  these  models  (Fel'zenbaum  1965,  Campbell  1965,  Coonetal  1974),  the 
compactness  N  of  pack  ice,  the  amount  of  ice  per  unit  area,  has  been  disregarded,  while  others 
(Doronln  1970,  Drogaitsev  1956,  Nikiforov  1957,  Rothrock  1970)  considered  the  compactness 
only  approximately.  Since  the  compactness  is  important  in  the  study  of  pack  ice  movement,  par¬ 
ticularly  in  the  marginal  region,  it  is  useful  to  derive  the  exact  equation  of  compactness  of  pack 
ice  and  to  show  the  main  effect  of  compactness  on  the  ice  drift.  We  are  going  to  derive  the 
equation  of  compactness  of  pack  ice  based  on  the  two-phase  flow  model  of  Pai-Li  (1977). 

2.  Pai-Li  model  of  Arctic  Pack  Ice 

The  main  concept  of  Pai-Li  model  is  the  two-phase  flow  approach  in  which  the  Arctic  pack 

ice  is  considered  as  a  mixture  of  solids  (ice  floes)  and  a  fluid  (sea  water).  Let  the  depth  of  ice 

be  H(x,y,  t)  which  is  in  general  a  function  of  x  and  y,  the  two  spatial  coordinates  on  tbe  surface 

of  the  ocean,  and  the  time  t.  The  maximum  thickness  of  ice,  Hq,  is  of  the  order  of  ten  meters. 

We  choose  a  typical  length  of  the  ocean  L  to  represent  the  dimension  of  the  surface  of  ocean  which 

we  are  interested  in  such  as  x  ~L  and  y  ~L.  The  representative  length  L  may  be  of  the  order 

of  1000  km.  Hence  L  is  much  larger  than  Hq.  We  may  divide  the  ice  packed  ocean  into  two 

layers.  In  the  upper  layer,  i.e. ,  z  <Hq,  where  z  is  the  coordinate  perpendicular  to  the  surface 

of  the  ocean  measured  from  the  surface  of  the  ioe  positive  downward,  we  have  a  mixture  of  ice 

floes  and  sea  water.  Since  L  »H  ,  the  toe  floes  may  be  considered  as  small  solid  particles  in 

o 

a  fluid  (sea  water).  In  this  mathematical  model,  the  concept  of  fluidization  of  solid  partloles  in 
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a  fluid  (Pal,  1971)  Is  used.  Hence  these  Ice  floes  may  be  considered  as  a  pseudo-fluid.  For 
z  4  Hq,  we  have  a  mixture  of  two  fluids:  one  Is  the  sea  water  and  the  othei  Is  the  pseudo-fluid 
of  ice  floes.  We  may  extend  the  theory  of  two-phase  flows  of  a  mixture  of  small  particles  and 
a  fluid  to  study  the  fundamental  equations  of  dynamics  of  pack  ice  (Pai-Li,  1976). 

Because  of  the  fact  that  L  »Hq,  we  may  reduce  the  three-dimensional  ice  dynamics  to 
an  effective  two-dimensional  problem  by  integrating  the  three-dimensional  equations  of  ice 
dynamics  with  respect  to  z,  the  vertical  coordinate,  with  the  limits  z  =  O  to  z  =  H(x,y,  t).  The 
resultant  equations  depend  only  on  two  spatial  coordinates  x  and  y  and  time  t  with  the  thickness 
of  ice  H(x,y,  t)  as  an  unknown  parameter  to  be  determined.  All  the  mathematical  models 
(Fel'zenbaum  1958,  Campbell  1965,  Doronln  1970,  Coon  et  al  1974)  are  based  essentially  on 
this  effective  two-dimensional  approach,  even  though  no  one  (except  Pai-Li  1976)  has  pre¬ 
viously  presented  the  fundamental  equations  of  the  effective  two-dimensional  ice  flow  explicitly 
in  this  manner. 

For  z  >H,  the  two-phase  flow  equation  will  be  reduced  to  those  of  ocean  water  when  the 
compactness  is  zero. 

3.  Derivation  of  the  Equation  of  Compactness  of  Pack  Ice  (Pai-Li  1976). 

The  equation  of  continuity  of  the  ice  of  the  three-dimensional  flow  in  the  two-phase  flow 
theory  of  ice  dynamics  is  as  follows: 
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where  Pl  =  izi  is  the  partial  density  of  the  pack  ice.  Pi  is  the  species  density  of  the  pack 
ice  and,  Z{,  is  the  volume  fraction  of  the  ice  in  the  mixture. 

For  two-phase  flow  of  a  mixture  of  ice  and  water,  we  have  two  different  definitions  of 
density  for  each  species  in  the  mixture.  Let  V  =  +  V  be  the  elementaiy  volume  of  the  mix¬ 

ture  where  is  the  volume  occupied  by  the  ice  and  is  the  volume  occupied  by  the  sea  water. 

Let  M  =  M,  +  M  be  the  total  mass  of  the  mixture  of  ice  and  water  in  V  where  M,  is  the  mass 
i  w  l 

of  ice  in  V  and  is  the  mass  of  water  in  V.  The  species  density  of  ice  is 


p,  -i 

v, 


=  Pt  (Tj*  p,  St,  c  ) 
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where  p{  is  in  general  a  function  of  ice  temperature  T^,  ice  pressure  P^,  ice  salinity  and 
ice  porosity  c  . 

The  partial  density  of  ice  P^  in  the  mixture  which  is  the  value  of  density  used  in  the  two- 
phase  flow  theory  (Pal,  1971),  i.e.,  the  Eq.  (1)  is 


1 


t 

1 

[ 

i 


i 

i. 

I 

to  • 

I 

I 

I 

i 

I 


2 


* 


p 


i 


where 
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(3) 


V 

Z.  =  1  =  volume  fraction  of  Ice  In  the  mixture.  (4) 

1  V 

Similarly,  we  have  species  density  of  sea  water  and  partial  density  of  sea  water  In  the 
mixture. 

The  source  term  of  Ice,  o^i  depends  on  the  mass  transfer  by  thermal  and  mechanical 
processes. 

Now  we  Integrate  Eq.  (1)  with  respect  to  z  from  z  =  O  to  z  =  H(x,y,  t)  and  have  the  follow¬ 
ing  equation 
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Since  the  ice  remains  on  the  top  layer  of  the  ocean,  we  may  assume  that  at  z  =  O  and  z  = 

H(x,y,  t)jWj  =  O,  or  if  we  consider  freezing  and  melting  of  ice,  the  rate  of  Increase  or  decrease 
at  the  surface  z  =  O  and  z  =  H  may  be  assumed  approximately  equal.  As  a  result,  the  fourth 
term  on  the  left-hand  side  of  Eq.  (5)  vanishes.  We  have  also  the  following  relations 


o  o 

and  similar  relations  for  the  second  and  the  third  terms  on  the  left-hand  side  of  Eq.  (5). 
Furthermore,  since  H(x,y,t)  is  much  smaller  with  respect  to  the  horizontal  dimension  L,  we 
may  assume  that  the  variation  of  p^,  Z{,  u^,  and  v^  with  respect  to  z  are  negligibly  small  and 
we  may  replace  their  value  by  average  values  of  these  variables  between  the  limits  z  =  O  and 
z  =  H  at  any  given  point  x  and  y  and  time  t.  Hence  we  have  the  following  relations 
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where  N  is  the  compactness  of  the  pack  ice  which  is  the  average  value  of  volume  fraction  of 
the  pack  ice  between  the  limits  z  =  Otoz=Hat  any  point  on  the  surface  of  the  ocean  (x,y) 
and  at  a  given  time  t.  All  the  average  variables  u^,  v^,  p^,  N  and  o{  (or  E^)  are  in  gen¬ 
eral  functions  of  x,  y  and  t  only. 

Substituting  Eqs.  (6)  and  (7)  in  (5),  we  have  the  equation  of  continuity  of  pack-ice  in 
polar  ocean  as  follows: 


3Np  jH  +  aNp^H  =  dNpjVjH 
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Equation  (8)  may  be  written  as  an  equation  for  compactness  N  as  follows 


It  is  interesting  to  notice  from  Eq.  (11)  that  the  compactness  equation  is  independent  of  the 
thickness  of  the  ice  H(x,y,  t).  Such  a  result  is  consistent  with  our  assumption  that  L  »H. 
This  result  is  essential  in  that  the  pack  ice  may  be  considered  as  a  two-dimensional  con¬ 
tinuum. 


4.  Comparison  of  our  Equation  of  Compactness  of  Pack  Ice  With  Similar  Equations 
in  the  Literature 

It  is  interesting  to  compare  our  equation  of  compactness  of  pack  ice  N,  (10)  or  (11) 
with  those  equations  of  compactness  in  the  literature. 

Rothrock  (1970)  gave  the  following  equation  of  compactness  [Eq.  (4)  of  Rochrock,  1970]: 


PtN  V_  =  -  N  p,Hft  +1V\  “  N  p .H  iji  +  *  (12) 

Dt  1  \dx  ay  /  1 

The  differences  between  (10)  and  (12)  are  tri-fold: 
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(i)  On  the  left-hand  sides  of  these  equations,  H  in  (10)  replaoes  N  of  Roth  rock's 

(12) .  It  Is  evident  that  N  In  (12)  of  Rothrook  must  be  a  printing  error  beoause  with  N,  die 
dimensions  of  (he  term  are  not  die  same  as  those  of  all  the  other  terms  In  (12). 

(11)  NH  Dj  Pt/Dt  In  Eq.  (10)  replaoes  N  p^H<|>  of  Eq.  (12)  of  Rothrook  who  had  little 

information  about  his  function  <|i .  We  have  no  difficulty  In  explaining  our  term  NHD{  pf/Dt. 

(ill)  There  are  some  basic  differences  between  Eqs.  (10)  and  (12)  In  concept.  In 

general,  the  velocities  of  ice  (u.  and  v.)  are  different  from  those  of  water  (u  and  v  )  In 

'll  '  w  w 

the  top  layer  O  H.  The  distinction  between  these  velocities  should  be  noted  and  In  Eq. 
(10),  the  use  of  operator  D{/Dt  should  be  kept  In  mind.  In  Rothrock's  analysis  (1970), 
such  a  distinction  had  not  been  mentioned.  Rothrook  separated  the  source  term  Into  two 
parts  and  in  (12),  only  one  part  of  the  source  term,  l.e. ,  the  rate  of  formation  of  ice  by 
freezing  open  water  is  included.  In  Rothrock's  analysis,  if  one  makes  the  distinction  be¬ 
tween  freezing  at  the  surface  z  =  O  and  freezing  at  O  <z  <H,  the  variation  with  z  of  all 
variables  should  be  considered  and  one  cannot  use  the  effective  two-dimensional  approach. 
For  effective  two-dimensional  approach,  the  source  term  should  be  the  average  total  source 
volving  the  mixture  of  ice  and  water  at  the  point  (x,y)  and  time  t.  Hence  ourE^  =  H  o{ 
is  different  from  Rothrock's  $  . 

A 

Various  Russian  authors  used  simplified  equations  for  compactness  N.  We  are 
going  to  compare  their  simplified  equations  with  our  exact  equation  of  compactness  (8), 

(10),  or  (11). 

Drogaltsev  (1956)  considered  the  conservation  of  the  mass  of  Ice  per  unit  area  and 
obtained  the  following  equation  of  compactness  for  two  dimensional  Ice  dynamics: 

*N  P,H  ♦  al<PlH"l  +«”P.HT1  =  0  ,13, 

dt  dx  dy 

Comparing  Eq.  (13)  with  our  exact  equation  (8),  we  see  that  Drogaltsev  neglected  the 
source  term  E^.  Furthermore,  he  missed  the  important  term  due  to  the  variation  of  the 
thickness  of  ice  H.  This  is  the  main  difference  between  the  effective  two-dimensional 
treatment  of  Ice  dynamics  and  the  ordinary  two-dimensional  approach  obtained  by  simply 
dropping  the  variation  with  z -coordinate.  For  lnstanoe,  If  we  drop  the  z-varlation  in  the 
three-dimensional  equations  of  motion,  we  will  never  get  the  Important  terms  due  to  air 
and  water  stresses  on  the  surfaoe  of  the  Ice.  Thus  such  a  simple  equation  of  continuity 

(13)  Is  not  consistent  with  the  correct  effective  two-dimensional  equations  of  motion  (15) 
and  (16).  Therefore  (13)  is  oorreot  only  If  we  neglect  the  source  term  and  if  we  assume 
that  Hie  thickness  of  ice  H  Is  a  constant 
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Nikiforov  (1957)  and  Doronln  (1970)  used  the  following  equation  for  compactness: 

_3N  =  -  /3NUi  +  J?Vl\  (14) 

at  \  ax  ay  / 

If  we  compare  Eq.  (14)  with  (11),  we  see  that  both  Nikiforov  and  Doronln  neglected  the  source  ( 

term  o{  and  the  variation  of  the  species  density  of  ioe.  Hence  (14)  is  valid  only  if  we  assume 
that  there  is  no  source  term  and  that  the  variation  of  species  density  of  ice  p^  is  negligible.  j 

5.  Effect  of  Compactness  of  Ice  on  Ice  Drift 

t 

In  studying  the  effect  of  compactness  of  ice  on  ice  drift,  we  should  solve  the  equation  of  i 

compactness  together  with  the  equations  of  motion  of  ice  and  the  equations  of  motion  of  air 

\ 

above  the  ice  layer  and  those  of  water  under  the  ice  as  shown  in  references  by  Fel’zenbaum  1 

(1958)  and  Campbell  (1965).  The  effective  two-dimensional  equation  of  motion  of  ice  according 
to  Pai-Li  model  (1976)  including  the  compactness  of  ice  are  as  follows:  1 . 
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where  f  =  2ft  sin$  =  coriolis  force  factor,  $  is  the  latitude,  and  fi  is  the  angular  velocity  of 

the  Earth.  The  terms  t,  and  t,  are  respectively  the  xz  and  yz  wind  stresses  on  the  top 

ixzo  iyzo 

of  the  surface  of  ice  floes,  z  =  O  and  the  terms  t,  „  and  Xt  are  respectively  the  xz  and  j 

ixzH  iyzH  £ . 

yz  water  stresses  on  the  lower  surface  of  the  ice  floes  z  =  H. 

There  are  two  types  of  interaction  forces  in  (15)  and  (16):  (i)  One  is  the  friction  force  due  * 

i  ■ 

to  the  difference  in  velocities  between  ice  and  water  i.  e. ,  slip  velocity  (u  -  u  )  and  the  other  is 

1  w  , 

the  internal  stresses  of  ice  floes  due  to  their  random  motion,  t.  ,  t,  ,  and  t  .  .  We  shall  ' 

ixx  ixy  iyy  { 

discuss  the  solution  of  Eqs.  (10),  (15)  and  (16)  in  another  paper.  Here  we  shall  derive  a  formula 
for  wind  ice  drift  which  shows  the  effect  of  compactness  for  a  simplified  case.  j 

Similar  to  the  solution  of  Fel’zenbaum  (1958),  we  neglect  the  Inertial  term  as  well  as  the 
interaction  force  terms.  Then  the  basic  equations  of  motion  are  reduced  to  a  set  of  linear  dif-  | 

ferentlal  equations  for  which  the  principle  of  superposition  is  applicable.  In  other  words,  as 
shown  by  Fel’zenbaum,  we  may  divide  the  ice  drift  into  two  parts:  One  is  the  pure  wind  drift  | 
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and  the  other  is  the  gradient  drift.  As  we  intend  to  compare  our  results  with  the  U.  S.  Naval 
Fleet  Weather  Facility  ice  drift  model,  and  their  model  is  for  wind  ice  drift  only,  we  shall  con¬ 
sider  only  the  pure  wind  drift  caused  directly  by  the  driving  effect  of  wind  on  ice  in  this  section. 
In  practice,  the  effect  of  permanent  currents  is  added  vectorially  to  the  wind  drift  to  determine 
the  total  drift  of  pack  ice. 

Similar  to  Fel' zenbaum' s  simplified  solution  for  the  modulus  of  the  pure  wind  drift  com¬ 
ponent  u  for  the  deep  water  part  of  the  Arctic  Ocean,  we  have  the  following  formula: 


(17) 


y  1  +  2m  +  2m 


where  po  =  density  of  air  ^  0. 0013  gr/cc;  p^  =  density  of  sea  water  as  1. 028  gr/cc;  A&  =  turbu= 


w 


lent  exchange  coefficient  of  air;  A  =  turbulent  exchange  coefficient  of  sea  water,  and  V  is  the 


w  G 

geostropic  wind.  Eq.  (17)  may  be  derived  from  Eq.  (58)  of  Fel'senbaum's  paper  (1958k  By 


comparing  our  equations  of  motion  of  ice  with  Fel* zenbaum' s  equations,  we  conclude  the  Eq.  (17) 
may  be  derived  from  our  equations  (15)  and  (16)  under  similar  approximation  provided  that  we 
replace  p ^  by  N  p^.  In  other  words  the  parameter  m  for  our  theory  should  be  defined  as  follows: 

_  Np, 


m  = 


ah 


(18) 


rw 


where  h  is  the  thickness  of  ice  and  a 


a  =  H  (19) 

V2AW 

It  ia  understood  that  in  Fel'zenbaum's  analysis,  N  =  1  is  implicitly  assumed. 

Even  though  for  actual  prediction  of  ice  drift  we  have  to  solve  the  fundamental  equations  of 
ice  dynamics,  we  may  estimate  some  first  order  effects  of  conpactness  by  considering  some 
typical  values  of  the  following  quantities: 

According  to  Fel'zenbaum  (1958),  we  may  assume 


and 


p{ah 


150 


1 


(20) 


(21) 
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Substituting  Eqs.  (18),  (20)  and  (21),  into  (17),  we  have 


u 


ia  =  0. 0154 


=  G(N).  V 


G 


(22) 


(1  +2N  +  2N 

when  G(N)  is  the  pure  wind  drift  factor  defined  in  Eq.  (22).  Equation  (22)  shows  some  essential 
effect  of  compactness  on  pure  wind  drift.  The  main  effect  of  compactness  is  that  as  compactness 
N  decreases,  the  pure  wind  drift  increases. 

It  is  interesting  to  compare  Eq.  (22)  with  the  empirical  formulas  of  pure  wind  drift  as  used 
currently  by  the  U.S.  Naval  Fleet  Weather  Facility  (FWF)  sea  ice  division  as  follows: 

The  empirical  formulas  used  by  FWF  are 


u.  =0.02  V  =  0.  02  v  for  C  2  knots 
la  s  g 


u.  =  0. 025  +  0. 00772  V  for  V  >  2  knots 
ia  G  G 


(23a) 

(23b) 


where  Vg  is  the  surface  wind  velocity. 

Equation  (23a)  is  essentially  Zubov's  formula  (1945).  The  "exact?'  coefficient  based  on 
Zubov’s  and  Nansen's  data  (Zubov  1945)  is 


u,  =  0.0155  V  =  0.0077  V 
xa  s  G 

But  Zubov  (1945)  in  his  book  rounded  the  figures  and  wrote 


(24) 


u,  =  0.02  V  =  0.01  V 
la  s  G 


(24a) 


where  Zubov  assumed  2V  =  V  .  We  should  compare  our  formula  (22)  with  the  "exact"  formula 

8  G  — ,  — ,  — » 

of  Zubov  (24).  FWF  assumed  that  for  low  speed,  V  <  2  knots,  V  =  V  .  Hence  they  have 

G  s  G 


u.  =  0. 0155  V  —  0. 0155  V  3  0. 02  V 
ia  s  G  G 


(24b) 


Equation  (23b)  is  based  on  the  empirical  formula  by  Skiles  (1968).  The  slope  of  the  line  u  vs 

-*  -»  -* 

VG  °*  *s  eclual  to  Zubov’s  value.  If  we  assume  that  the  slip  velocity  (u^  -  u^)  in  Eqs.  (15) 

and  (16)  as  a  constant,  we  may  obtain  the  constant  term  in  (23b)  by  choosing  proper  friction  co¬ 
efficient.  We  shall  discuss  this  point  in  another  paper. 

From  (22)  we  have  the  following  values: 

N  =  0. 1  0.3  0.5  0.7  0.8  1 

G(N)  =  0.0139  0.0115  0,0097  0.0084  0.0077  0.0069 

Our  formula  (22)  is  identical  to  Zubov's  formula  (24)  if  we  take  N  3  0.  8. 

Pure  wind  drift  would  be  deflected  from  the  isobaric  drift  to  the  right  at  an  angle  P  given  by 
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Fel'zenbaum  [(59)  of  Fen'zenbaum  1958] 

tan  p  =  1  -  1  (25) 


pjNha 


Again,  In  our  theory,  we  may  replace  p ^  by  Np^  In  Fel'zenbaum' s  formula.  By  the  approxima¬ 
tion  (21),  (25)  becomes 

tan  P  =  1  (26) 

1  +1 
N 

Equation  (26)  gives  the  following  values 

N  =  0. 1  0.3  0.5  0.7  0.8  1 

P  =  5  “12*  13°  18°26*  22*22*  24°18»  26*34* 

Both  (22)  and  (26)  show  qualitatively  the  effect  of  compactness  N  on  the  pure  wind  drift 
and  the  deflection  angle  of  pure  wind  drift  from  lsobasic  drift.  However,  because  the  approxi¬ 
mation  of  (21),  we  do  not  expect  that  these  formulas  give  exact  quantitative  results  from  field 
observation.  For  using  these  formula  to  check  the  field  observations,  we  should  take  the  actual 
variation  of  ice  thickness  Into  consideration. 

6.  Conclusions  and  Recommendations 

(1)  We  have  given  the  accurate  equation  of  compactness  of  pack  ice  in  Arctic  Ocean  based 
on  Pai-Li  model  (1976)  of  two-phase  flow  theory  and  compare  it  with  the  simplified  equations  of 
compactness  of  pack  ice  in  literature. 

(2)  A  simple  formula  for  pure  wind  ice  drift  with  compactness  effect  is  derived.  This 
formula  agrees  well  with  empirical  formula  of  Zubov  and  those  used  by  U.S.  Naval  Fleet  Weather 
Facility  for  ice  drift  prediction.  Our  formulas  show  qualitatively  the  effect  of  compactness  on 
pure  wind  ice  drift  and  the  deflection  of  pure  wind  ice  drift  from  lsobarlc  direction.  As  the  com¬ 
pactness  decreases,  the  pure  wind  ice  drift  increases  for  a  given  geostropic  wind  and  the  de¬ 
flection  angle  decreases. 

(3)  We  shall  use  our  fundamental  equations  of  ice  dynamics  based  on  Pai-Ll  model  of  two- 
phase  flow  theory  for  the  calculation  of  ice  drift  in  Arctic  Seas,  particularly  in  the  marginal 
sea  area  where  the  effect  of  compactness  of  ice  would  be  Important. 
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